Morphology of the Superconducting Vortex Lattice in Ultra-Pure Niobium 
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The morphology of the superconducting flux line lattice (FLL) of Nb comprises gradual variations 
with various lock-in transitions and symmetry breaking rotations. We report a comprehensive small- 
angle neutron scattering study of the FLL in an ultra-pure single crystal of Nb as a function of the 
orientation of the applied magnetic field. We attribute the general morphology of the FLL and its 
orientation to three dominant mechanisms. First, non-local contributions, second, the transition 
between open and closed Fermi surface sheets and, third, the intermediate mixed state between the 
Meissner and the Shubnikov phase. 

PACS numbers: 71.18.+y, 74.25.Dw, 74.25. Q+, 78.70.Nx 



In the conventional Ginzburg-Landau theory of super- 
conductivity the formation of a flux line lattice (FLL) 
- the so-called Shubnikov phase - is the result of repul- 
sive flux line interactions, where the coherence length is 
a local property of the superconducting state. Numerous 
experimental studies have shown [lEl |3] , that the FLL 
depends sensitively on the microscopic nature, symme- 
tries and number of superconducting gaps, as well as the 
general topology of the Fermi surface and the effects of 
impurities. This raises the question how to generalize the 
description of FLLs and makes studies of conventional 
superconductors of general interest. 

The superconductivity in Nb is ideally suited to pro- 
vide such general information ^ i5j. The properties of 
Nb are characteristic of a single gap [6l [7] , avoiding the 
complexities of multi-gap systems [U [H |9] . The low value 
of the Ginzburg-Landau parameter, n ^ 0.8, places Nb 
at the immediate border between type 1 and type 2 su- 
perconductivity. This implies that the lower critical field 
jiioHci{T ^ 0) = 120 mT is high and the intermediate 
mixed state (IMS) between the Meissner and the Shub- 
nikov phase is extended, thus making the IMS and the 
underlying change from attractive to repulsive vortex in- 
teractions experimentally accessible [10]. Moreover, the 
large coherence length of Nb implies that the Fermi sur- 
face topology and thus non-local effects are important. 

For magnetic field H parallel to a crystallographic 
(111) direction the FLL in Nb exhibits a hexagonal sym- 
metry. In contrast, for H \\ (100), a rich phase di- 
agram unfolds shown in Fig.jljA), which prior to our 
study was understood as follows [H O [11]: A two- fold 
isosceles phase and three distinct FLL phases are ob- 
served in the vicinity of Tc and at low temperatures, re- 
spectively [H O HI]. For increasing H at low-T these 
are at first a square phase (low- field square), a scalene 



phase at intermediate fields and a square phase near Hc2 
(high- field-square). Based on the magnetic field depen- 
dence for (100) it had been speculated that the low- field- 
square may be related to the IMS [4 . It had further 
been established that the three low-T phases are tilted 
with respect to the four-fold symmetry of the underly- 
ing crystal structure into a low symmetry direction, thus 
additionally breaking the crystal mirror symmetry. The 
two-fold isosceles phase in comparison breaks crystal ro- 
tational symmetry but still obeys mirror symmetry. It 
had finally also been established that the FLL symme- 
tries change and the tilting vanishes when rotating the 
magnetic field away from (100) until the equilateral FLL 
is recovered for (111). 

However, in these studies the precise evolution of the 
FLL morphology as a function of field orientation, nec- 
essary to identify the nature of the morphology and tilt- 
ing, remained open. This originated in incomplete data 
sets and, more importantly, in the purity of the samples 
studied so far, which still displayed distinct pinning ef- 
fects despite their residual resistivity ratios RRR ^ 1000. 
In this Letter we report a systematic study of the FLL 
symmetries using small angle neutron scattering (SANS). 
To resolve the open issues we have studied an ultra-pure 
single crystal of unprecedented purity with RRR > 10^, 
which displays essentially no pinning. 

Our study identifies the low-field-square as a property 
of the IMS. The remaining FLL phases may be viewed as 
driven by a gradual evolution from an equilateral FLL at 
high temperature and low field to a square FLL at low 
temperature and high field (cf Fig.[l](B)), with a lock- 
in transition below a magic angle near (100), where the 
locked phases are also tilted. The scalene FLL thereby 
appears in transitional regions of high frustration. Our 
experimental results show strong analogies with an evolu- 
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FIG. 1: (A) Schematic flux lattice phase diagram for a nio- 
bium sample with a finite demagnetizing factor in reciprocal 
space with the magnetic field applied along a four-fold (100) 
axis, where 77 = 0° denotes the angle between the applied 
magnetic field and the (100) axis in a (Oil) plane. (B) Phase 
diagram for 77 = 25°. See text for further details. 

tion of non-local corrections in an Eilenberger [12 treat- 
ment of a Fermi surface with four-fold symmetry [13] for 
high-/^ superconductors. The magic angle of the transi- 
tion to tilting corresponds remarkably well to a change 
between open and closed Fermi surface sheets, empiri- 
cally suspected to affect the superconductivity [14 . 

For our studies, a cylindrical Nb single crystal with 
a length of 20 mm and a diameter of 4.5 mm was cut 
from a rod that had been produced at the ZFW Dres- 
den [15 more than 30 years ago. The preparation pro- 
cess consisted of purification by liquid-liquid extraction 
combined with chlorination and thermal decomposition 
of NbCls followed by electron beam floating zone melt- 
ing, decarburization in oxygen atmosphere and annealing 
in UHV [m [16] . The impurity content was estimated to 
be less than 1 ppm for interstitial and better than 2 ppm 
for substitutional impurities. The sample was again RF- 
annealed in UHV above 2000° at the University of Birm- 
ingham for one week to remove interstitials followed by 
surface oxygenation to reduce the Bean-Livingston bar- 
rier for surface pinning [4 . The RRR was measured with 
an eddy current decay method at the University of Birm- 
ingham and found to exceed ~ 10^. The AC susceptibil- 
ity and the magnetization, measured at the Technische 
Universitat Miinchen were consistent with the literature. 

The measurements were carried out at the SANS 
diffractometer V4 at BENSC. Neutrons with a wave- 
length A = 12 A ±5.5% were collimated over a distance 
of 12 m before reaching the sample, with a multidetec- 
tor at a distance of 12 m. An 'orange' cryostat was used 
for measurements down to 1.5 K and a horizontal magnet 
for fields up to 500 mT. The sample was mounted on a 
motorized sample stick that could be rotated about the 
vertical axis with an accuracy of ±0.05° (cf. Fig.|2|A)). 
The cylinder axis of the sample coincided with the rota- 
tion axis and corresponded to a crystalline [Oil] direc- 
tion. Magnetic fields were applied in the (Oil) plane, 
perpendicular to the cylinder axis so that demagnetizing 
effects did not change during our measurements. Thus 
the direction of H could be rotated away from the [100] 
axis towards a [111] axis in the (Oil) plane [17]. Akin 
phenomena to those reported here have been obtained 
for a (100) plane of rotation on a different sample [5l. In 




FIG. 2: (A) Small angle neutron scattering set-up. (j) denotes 
the angle between the incoming neutron beam and the mag- 
netic field, while 77 denotes the angle between magnetic field 
and the crystallographic (100) direction. (B) Typical detector 
image in the scalene FLL phase. For clarity just one domain 
is marked in red. (C) The nomenclature for describing the 
FLL: a, f3 and 7 represent the internal angles of the FLL unit 
cell, S represents the tilting angle of the FLL with respect to 
the horizontal [111] axis. The real space image of the FLL 
can be obtained by a rotation of the scattering patterns (B) 
and (C) by 90° around the neutron beam axis. 

the following we denote with (j) the rocking-angle between 
the incoming neutron beam and while rj denotes the 
angle between H and the (100) direction. 

For the precise sample alignment within ±0.2° we ex- 
ploited the degeneracy of the superconducting domains 
for field parallel (100) as described in Ref. [5^. Rocking 
scans of (j) were performed for each measurement point. 
To avoid any hysteretic effects, each data point was taken 
after field cooling the sample [17 . Several domain popu- 
lations may be seen for 77 = at jhqH = 250 mT and 
T = 1.7 K (Fig.[2|B)), as expected for the symmetry 
breaking of the FLL with respect to the four-fold sym- 
metry of the crystal lattice. For clarity only one of these 
domains has been marked with red lines. To extract the 
morphology of the FLL the scattered intensities were fit- 
ted by 2D-Gaussians. As shown in Fig.|2|B)(C) a, P and 
7 describe the internal angles of the scattering pattern, 
defining the structure of the primitive unit cell, while S 
describes the orientation of the primitive FLL unit cell 
with respect to the horizontal [111] axis. Previous stud- 
ies have shown that changes of domain population of the 
FLL already occur within small changes of 77, while the 
shape and orientation of the unit cell are conserved 0. 
Thus we focus on the morphology of the FLL. 

The field dependence of the integrated scattering in- 
tensity of a FLL Bragg spot and its reciprocal lattice 
spacing \Q\ are shown in Fig. [3] for T = 3.6 K and 7^ = 0°. 
We denote the constant flux density due to the attractive 
flux line interactions in the IMS with Hq. As expected 
for our cylindrical sample with demagnetization factor of 
^, no intensity is observed for both increasing and de- 
creasing magnetic field up to ^/iqHci = 65 mT, followed 
by a gradual increase up to ^iiqHq ± ^/j^qHci ^ 115 mT 
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FIG. 3: Integrated scattering intensity over a rocking curve 
and reciprocal lattice vector of the SANS diffraction pattern 
as a function of increasing and decreasing applied magnetic 
field in the isosceles FLL phase. Only one of two characteristic 
values of IQI is shown, as only rocking curves about a vertical 
axis could be carried out. Note the complete Meissner flux 
expulsion for both increasing and decreasing field. 




FIG. 4: Evolution of the flux lattice as a function of 77 for 
various temperatures and magnetic fields, where typical sums 
over rocking scans are shown for the high- field-square (i), sca- 
lene (ii), low field square (iii) and isosceles (iv) structures. All 
FLL structures are given in reciprocal space, with the angles 
plotted in the l.h.s. panels defined in Fig[2] The scattered 
intensity is plotted on a logarithmic scale. Due to a slight 
misalignment of the sample, the scattering patterns are ro- 
tated by ~ 1° around the beam axis 17 . 



[18]. To the best of our knowledge such a complete 
Meissner fiux expulsion for decreasing H has never been 
seen before in SANS, underscoring the outstanding purity 
and perfection of our sample. Above ^iiqHq + ^ji^Hd 
the intensity decreases characteristic of a conventional 
Shubnikov phase. The value of Q reversibly assumes 
a constant value Q ^ 0.004 between ^/io^ci and 
^fioHo -\- ^fioHci, identifying this regime as the IMS. In 
contrast the integrated intensity displays hysteretic be- 
havior between 110 mT and 125 mT. Moreover, also the 
rocking width of the FLL is increased for decreasing field. 
The same behavior in sign and magnitude has been ob- 
served in the IMS in low quality Nb [18] . The data in our 
high quality sample hence reflect an intrinsic property, 
most likely due to dendritic flux intrusion for increasing 
fields in contrast to FLL islands for decreasing fields [19] . 

We have studied the evolution of the FLL morphology 
in detail for 1.5K < T < 5.5 K, 110 mT < H < 330 mT 
and ?7 < |±25°|. This covered the parameter range nec- 
essary to determine the FLL phase diagram. Shown 
in Fig.[4](i) are typical data of the high-field-square for 
77 = 0° with a = /3 ^ 45°, 7 ~ 90° and a tilting angle 
S ^ 11° . As a function of rotation angle rj the symmetry 
breaking vanishes for rjc ^ 15°. Above rjc the internal 
angles a = P increase and 7 decreases characteristic of 
an isosceles lattice that transforms into the hexagonal 
lattice with a = /3 = 7 = 60° for H || (111). 

Fig.[4](ii) shows typical data in the scalene phase, 
which is characterized hy a ^ 50° , P ^ 60° and 7 ^ 70° . 
Like the high-field-square the scalene phase is tilted, 
where S ~ 14°. The tilting vanishes at rjc ^ 15°, where 
the scalene structure turns isosceles, i.e., a = /3 ^ 50° 
and 7 ~ 80°. Approaching (111) the FLL turns hexag- 
onal. Typical data of the low- field square are shown in 
Fig.[4](iii) (a = /3 = 45°, 7 = 90°), which abruptly turns 
into a scalene lattice already for an angle rji 5°, while 
the tilting of the flux lattice exists up to r]c ^ 20°, i.e., in 
contrast to all other phases the fundamental symmetry 
and the tilting exhibit distinctly different critical angles. 

Finally, the reorientation angle rjc ^ 20° is also im- 
portant in the isosceles phase at high temperatures {a = 

^ 62°) as shown in Fig.|4](iv) [U]. For 7^ = 0° the 
unit cell is orientated such that a [010] axis acts as mir- 
ror plane, yielding f + = 45°, i.e., 13°. The 
two-fold isosceles phase hence does break the underly- 
ing four-fold rotational symmetry of the lattice at 77 = 0° 
but still obeys mirror symmetry. With increasing rj this 
orientation changes between rj ^ 10° and rjc such that 
[Oil] becomes a mirror plane of the isosceles unit cell 
(J=0). During this reorientation the unit cell is scalene. 
This may be compared with the low-field-square shown 
in Fig.|4](iii), which becomes scalene already at rji ^ 5° 
as discussed above. At rjc ^ 20° both the tilting vanishes 
and the lattice returns to the isosceles structure. 

Thus rjc represents a magic angle that is essentially 
the same for all FLL phases (Fig.[5](C) and (D)). This is 
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FIG. 5: Evolution of the internal angles a,P and 7 for 77 = 20° 
as a function of magnetic field for low temperatures (A) and 
as a function of temperature for 100 mT and 150 mT (B). Note 
the trend towards a square lattice for low T and high H. Pan- 
els (C) and (D) show the temperature and field dependence 
of ?7c, which is essentially unchanged, where the error bars 
indicate the range of the transition. 

strongly contrasted by r]i of the low- field-square. More- 
over, the reversible transition fields of the low- field- 
square and the IMS shown in Figs. [3] and [4] places the 
low- field-square right in the IMS. Due to the large ^gl 
as compared to the spacing of the flux lines, the interac- 
tion between flux lines in the IMS includes contributions 
from an overlap of the vortices. This leads to a long-range 
attractive FLL interaction in the IMS, where a square 
FLL is expected in a theoretical analysis, ignoring any 
symmetries imposed by the crystal lattice consis- 
tent with our experiment. Taken together this identifies 
the low-field square as a property of the IMS. 

With the low-field-square being a property of the IMS 
we may now address the remaining phase diagram. We 
first consider gradual variations of the FLL morphologies, 
where it is instructive to begin with the isosceles phase 
for T] > r]c (cf Fig.jljB)). As shown in Fig.jsjfor r] ^ 20° 
we find that a = /3 ^ 45° and 7 90° for decreas- 
ing temperature and increasing magnetic field. In other 
words the lower the temperature and the higher the mag- 
netic field the stronger the trend to form a square lattice. 
This is underscored by a weak four-fold anisotropy of Hd 
[TTJ [21] . The same trend to form a square lattice is also 
present for r] < r]c^ where the FLL morphology locks into 
distinct FLL phases. A scalene FLL thereby emerges al- 
ways in the transition regime between structures with 
higher symmetry (isosceles or square), i.e., under condi- 
tions of maximum frustration. This may be readily seen 
in Figs.[4](ii), (iii) and (iv), where we find that a and p 
differ for r] < rjc^ while a = /3 for r] > rjc. Likewise, a 
scalene FLL due to frustration is also seen during the 
rotation between the two isosceles phases in Fig.[4](iv) 
with [Oil] and [010] acting as mirror planes, respectively. 
Hence the same systematic trend is observed for all re- 
orientation transitions shown in Fig.|4](ii), (iii) and (iv), 
but where the width of the transition differs. 

In a pioneering study Nakai et al [13 have shown for 



high-/^ superconductors, that the trend observed here is 
expected for a four-fold Fermi surface symmetry and re- 
lated four-fold gap anisotropy. This reproduces qualita- 
tively the behaviour with H \\ (100), where a decreasing 
amount of this trend may be expected to survive when 
turning towards (111). 

We finally note, that roughly the same value of r]c and 
the size of the tilting, 6c, are observed across the entire 
phase diagram. These values are remarkably close to the 
angle separating open and closed Fermi surface sheets 
inferred from the magnetoresistance and Fermi surface 
calculations [22j |23] , empirically suspected to be impor- 
tant for the superconductivity [14 . Moreover, the tilting 
angle 5c compares very well with the size of certain necks 
in the Fermi surface [24 . Thus the trend of the FLL mor- 
phology to lock into stable phases wins at the expense of 
the symmetry breaking rotations as driven by the Fermi 
surface topology. 

In conclusion, we attribute the general morphology of 
the FLL in Nb and its orientation to three dominant 
mechanisms. First, non-local contributions as described 
by Nakai et al in the Eilenberger formalism [13 , second, 
the transition between open and closed Fermi surface 
sheets and, third, the intermediate mixed state at the 
border between the Meissner and the Shubnikov phase. 
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